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<N '. Abstract 
> 

10 I Nonequilibrium effects in layered superconductors forming a stack 

of intrinsic Josephson junctions are investigated. We discuss two basic 
00 . nonequilibrium effects caused by charge fluctuations on the supercon- 

I ducting layers: a) the shift of the chemical potential of the condensate 

and b) charge imbalance of quasi-particles, and study their influence 
on IV-curves and the position of Shapiro steps. 
I Pacs: 74.72.-h, 74.50.+r, 74.40. +k. Keywords: layered supercon- 

ductors, Josephson effect, nonequilibrium superconductivity 

■ 1 Introduction 

For the strongly anisotropic cuprate superconductors Bi2Sr2CaCu2 08+5 
^ ■ (BSCCO) and Tl2Ba2Ca2Cu30io+5 (TBCCO) the electromagnetic transport 

■ perpendicular to the CUO2 layers can well be described by a model where 

the Cu02 planes with the intermediate material form a stack of Josephson 
junctions. This intrinsic Josephson effect can be seen in the multibranch 
structure of the IV-curves but also in the behavior of the material in external 
magnetic fields and under high frequency irradiation |jl|, ||, 0, ^. Here also 
Shapiro steps have been observed ^ . 

In the presence of a bias current each junction of the stack can be either 
in the resistive or superconducting state. In the resisistive state a finite 
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Figure 1: Typical mesa-structure used for 2-point measurements consisting 
of a stack of superconducting layers n = 1,2..., a normal electrode n = on 
top and a base electrode. 



dc-voltage V appears together with voltage oscillations with a frequency uj 
given by the Josephson relation huj = 2eV. Both effects are accompanied 
with static and oscillating charge fluctuations on the layers. In a system of 
weakly coupled very thin atomic layers such charge fluctuations may lead 
to nonequilibrium effects described by a shift of the chemical potential of 
the condensate and a charge imbalance between electron and hole-like quasi- 
particles P, |^, ||, ||, |T0|, |lT], |l2l . This is different in classical Josephson systems 
with massive superconducting electrodes. Here similar charge fluctuations 
occur, however, they are restricted to the surface of the superconducting 
electrode and do not affect the superconductor as a whole and in particular 
do not reach the normal electrode. 

It is a non-trivial question wether such nonequilibrium effects have mea- 
surable consequences. In an experiment like that of Clarke |13| several con- 
tacts are necessary to measure the difference between the chemical potential 
of condensate and quasi-particles. Experiments on cuprate superconductors 
are mostly done on mesa structures (see figure 1), where the top layer is 
covered with one gold electrode, and the same electrode is used to inject the 
current and measure the voltage (for a 4-probe measurement see f^). Also 
in such situations nonequilibrium effects should be observable for instance 
through their influence on the IV-curves: there should be a difference in the 
total voltage depending on wether two resistive junctions are neighbors in the 
stack or separated by a junction in the superconducting state. Irregularities 
in the brushlike structure of the IV-curves, however, can also be produced 
by other effects, i.g. fluctuations in the critical current density. Therefore in 
this paper we investigate in particular the possible influence of nonequilib- 
rium effects on the voltage-position of Shapiro steps, which should be more 
precise. We also add some theoretical considerations for measurements on 
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systems with multiple contacts, where it is possible to measure the charge 
imbalance relaxation directly. 

In a Shapiro-type experiment high-frequency radiation is applied to the 
sample, in this case to the antenna connected to the gold contact on top of 
the mesa. For sufficiently strong radiation the internal Josephson oscillation 
frequency of a junction in the resistive state is locked to the external fre- 
quency uj fixing the voltage in a small current interval. For a single junction 
in the resistive state one expects a set of Shapiro steps in the IV-curve at 
voltages V = V^ont + mhuj /{2e). Here Kont is the voltage due to the contact 
resistance between the normal electrode and the first superconducting layer. 
It is determined by measuring separately the IV-curve in the superconduct- 
ing state. A similar result is expected, if m junctions are in the resistive 
state and we restrict ourselves to the first Shapiro step in each branch. In 
this paper we will show that the position of steps may be changed into 
y = Kont + mfujj/{2e) — 5V , where again V^ont is the contact resistance in 
the superconducting state. This happens if the first resistive junction is close 
to the normal electrode. The apparent reduction of the Shapiro step volt- 
age is indirectly due to a change of the contact voltage in the presence of 
a nonequilibrium charge distribution on the first superconducting layer, and 
5V is proportional to the relaxation time of the quasi-particle charge. A 
similar result is obtained for the position of the m*^ Shapiro step on the first 
branch. 

In a recent experiment on a BSCCO sample a downshift of the first 



Shapiro step on the first branch has been observed |]I4[. Experiments on 
other samples show more complicated structures which could not be inter- 
preted completely. Therefore in this paper we concentrate on the theoretical 
aspects of nonequilibrium effects in layered superconductors. In the follow- 
ing section we summarize the basic results of a theory developed by one of 
us (D.A.R.) [|ll|] which is based on a kinetic equation for the distribution 
function describing charge imbalance. First this theory is used to study the 
influence of nonequilibrium effects on the IV-curves. Then these results are 
applied to the calculation of Shapiro steps and flnally to the analysis of sys- 
tems with multiple contacts. In all these applications of the nonequilibrium 
theory we restrict ourselves to stationary effects. We also neglect the influ- 
ence of high frequency oscillations on the dc-component of the supercurrent, 
which is a reasonable approximation for weakly coupled superconducting 
layers (large values of the McCumber parameter). Some results from the 
numerical simulation of the time-dependent kinetic equations can be found 
TTl. 



m 



Let us briefly mention here also the relation of the present theory to other 
theoretical approaches. As mentioned above there are two nonequilibrium ef- 
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fects connected with charge fluctuations a) the shift of the chemical potential 
of the condensate, b) the charge imbalance of quasi-particles (there is also an 
influence of nonequilibrium effects on the amplitude of the superconducting 
order parameter, which will not be considered here). If we take into account 
the shift of the chemical potential but neglect charge imbalance we essen- 
tially obtain the results of the theory developed by Koyama and Tachiki 0. 
In the approach by Artemenko and Kobelkov and by our group ||10|, Q 
a systematic perturbation theory in the gauge-invariant scalar potential is 
used. In these theories charge imbalance is considered only indirectly as far 
as it is induced by fluctuations of the gauge invariant scalar potential. The 
present theory is more general. Here charge imbalance is taken into account 
as an indepent degree of freedom and therefore the results are different from 
earlier treatments. A detailed comparison of the different approaches will be 
given in a future publication. 



2 Theory of stationary nonequilibrium effects 
in layered superconductors 

2.1 Generalized Josephson relation 

We consider a system of superconducting layers n with superconducting order 
parameter A„(t) = |A| exp(zx„(t)) with time- dependent phase Xnif) neglect- 
ing a possible time dependence of the amplitude of the order parameter. The 
basic quantity which enters all the Josephson relations is the gauge invariant 
phase difference between layer n and n + 1 given by 

2e 

7n,n+l(t) = Xnif) - Xn+lif) - ^ J dzA^{z,t), (1) 

where Az{z,t) is the vector potential in the barrier. For the time derivative 
of 7n,n+i we obtain the generahzed Josephson relation: 

2e / \ 

7n,n+l = -J- [Vn,n+1 + $n+l " $nj ■ (2) 

Here 

rn+l 

Vn,n+1 = / dzEz{z,t) (3) 



is the voltage and ^n{t) is the so-called gauge invariant scalar potential de- 
fined by 

$n(t)=0n(t)-^Xn(t), (4) 
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where 0„(t) is the electrical scalar potential. In this paper e denotes the 
elementary charge. The charge of the electron is — e. The quantity /i7„,„+i is 
the total energy to transfer a Cooper pair between neighboring layers n and 
n + 1. The quantity e<l>„ can be considered as shift of the chemical potential 
of the superconducting condensate with respect to an avarage chemical po- 
tential. It determines the minimum in k-space of the quasi-particle excitation 
energy. 

In the following we have to consider fluctuations of the charge densi- 
ties on the superconducting layers. For this purpose let us briefly mention 
some basic concepts for the description of nonequilbrium superconductors 
as introduced by Tinkham and Clarke [0. In the BCS theory one may 
split the total electronic charge density on each layer into p = p'^ + with 

= J2k '^k P'^ = Y^i J2k{'^k~'^k)fk- In equilibrium the latter vanishes 
because there are equal numbers of electron-like and hole-like excitations. In 
a nonequilibrium state of the superconductor a finite quasi-particle charge 
may exist, which is called charge imbalance. The quantity p*, sometimes 
called condensate charge, should not be confused with the superfluid charge 
density, which vanishes at the transition temperature. On the contrary 
approaches the total density in the normal state. Fluctuations of the two 
charge densities are caused not only by shifts in the chemical potential en- 
tering the quasi-particle excitation energy and the functions Uk,Vk, but also 
by fluctuations in the distribution functions fk- In a quantum-kinetic theory 
based on quasi-classical Green's functions the same distinction between the 
two types of charge densities can be made. The k-dependent functions are 
then replaced by corresponding energy-dependent quantities. 

A fluctuation in the (2-dimensional) charge density p^ is directly related 
to the shift of the chemical potential in the layer n: 



Here A^(0) is the 2-dimensional density of states of the conduction electrons 
in the layer. Unlike a true chemical potential, which is deflned only in equi- 
librium, the quantity $n(^) ^lay be time dependent. In the following we will, 
however, restrict ourselves to stationary processes where $„ is independent 
of time. 

It is convenient to express also fluctuations in the charge imbalance with 
help of a quasi-particle potential by writing formally 



Pnit) 



2e^N{0)<l>n{t). 



(5) 



pl = 2e2iV(0)vl>„(t). 



(6) 




(7) 
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The calculation of the charge imbalance will be discussed in section (2.3) 
for the stationary case. 

With help of (J^) and the Maxwell equation [d is the distance between the 
layers) 

Pn = -J-iVn^n+l — Ki-l,n) (8) 

the generalized Josephson relation now reads: 

h 

— 7n,n+l = (1 + 2a)K,n+l " a(K-l,n + K+l,n+2) + " (9) 

with a = eQ/{2e^N{0)d). It shows that the Josephson oscillation frequency 
is determined not only by the voltage in the same junction but also by the 
voltages in neighboring junctions. Furthermore it is influenced by the quasi- 
particle potential \1/ on the layers. If we neglect the latter we obtain for 7n,n+i 
the same result as in 0. 

For a Josephson junction in the presence of sufficiently strong high- 
frequency irradiation the average Josephson oscillation frequency is locked 
to the external frequency u. From the generalized Josephson relation (P) we 
obtain for the m*^ Shapiro step 

= ^(7n,n+l) (10) 



;i + 2a)K,n+l - a{Vn-l,n + K+l,n+2) + - ^ 



11} 



where the dc-components of the quantities on the r.h.s have to be used. 
This relation can also be used to describe the voltage for a barrier in the 
superconducting state, then the l.h.s. is zero: 

= (1 + 2a)K,n+l - «(K-l,n + K+l,n+2) + ^„+l - ^n- (H) 



2.2 Current equation 

In order to calculate the electronic transport in layered superconductors we 
need expressions for the current and charge density. The results are derived 
with help of a Keldish technique for nonequilibrium Green's functions. In 
particular the quasi-particle charge is calculated from a kinetic equation for 
the corresponding quasi-particle distribution function. For the current den- 
sity the following approximate result is obtained: 

jn,n+l = 3c sin 7„,„+l + ^"'J"^^ (^7n,n+l + " ^n+l) • (12) 

The first term is the current density of Cooper pairs. The rest is the quasi- 
particle current density (we neglect here the interference term proportional 
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to COS 7) . The physical meaning of the latter becomes more apparent if we 
express these terms with help of (1^^) by the charge densities: 

which can also be written as: 

jn,n+i = jc sin 7„^„+i H '— ,n+l - a{Vn~l,n + Vn+l.n+2)) ■ (14) 

Then we see that the quasi-particle current is driven not only by the voltage 
but also by a diffusion term proportional to the charge difference in the two 
layers. The expression (|13D for the current density is valid for temperatures 
close to Tc. For T <^ the conductivities are different for the two contri- 
butions of the quasi-particle current (they correspond to the quantities ctq 
and (Ji introduced by Artemenko et al. §]). For simplicity we neglect this 
difference in the following. 



Note that ( p!3i |1J) stay meaningful also in the normal state for a system 
of junctions with different conductivities. Then a describes the shift of the 
chemical potential due to charge accumulation in the normal state. 

In the stationary state the current density jn,n+i is the same for all barriers 
and can be replaced by the bias current density j. If the Josephson junction 
between layers n and n + 1 is in the resistive state (without high-frequency 
irradiation), then we may use the current relation (|1^ with a vanishing dc- 
component of the supercurrent (which is a good approximation in the limit 
of large values of the McCumber parameter) and obtain: 

jd 

(1 + 2a)Vn,n+l - a(K-l,n + Vn+l,n+2)- (15) 



For a correct determination of the total voltage of a stack as shown in fig- 
ure 1 an investigation of the influence of nonequilibrium effects on the contact 
voltage between the normal electrode and the first superconducting layer is 
necessary. The physical nature of this contact is not known precisely. It may 
vary between a tunneling contact and a metallic contact with superconduc- 
tivity partially suppressed on the first superconducting layer by proximity 
effect. In the following we assume that this contact can be described by 
a tunneling contact between a thick normal electrode and a thin supercon- 
ducting layer. Then we may neglect the shift of the chemical potential in the 
normal layer and obtain for the current density between the normal electrode 
(n = 0) and the first superconducting layer {n = 1) (compare with (|l^)) 

[l + a)Vo,i-aV,,2. (16) 
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There might be a small shift of the chemical potential on a thin layer at 
the surface of the normal electrode. This can be included in the theory, 
but has no influence on the total voltage. For notational simplicity we have 
assumed here that the barrier width d between the normal electrode and the 
first superconducting layer is the same as that between two superconducting 
layers. This can easily be generalized but does not change the content of the 
final results. Note that the quasi-particle potentials appear explicitly only 
in the equations derived from the generalized Josephson equation and not in 
the current equations (|l^, p!6|). 

2.3 Relaxation of the quasi-particle charge 

In addition to the current equation the microscopic theory provides us with 
an equation for the charge density. As the quasi-particle charge is only part of 
the total density an additional mechanism for the charge imbalance relaxation 
has to be introduced. In the stationary case the following result for the quasi- 



particle charge is obtained ||11|| : 

p1 = -2e^A^(0)^„ = Tq{jc sin -fn,n+i -jc sin 7„_i,„), (17) 

where Tq is the charge imbalance relaxation time. 

This result can be made plausible in a simple way: Let us start from a 
balance equation for the superfiuid charge density: 

dp" 

= jcsin7„_i,„ - jcsin7„,„+i - -^\conv (18) 

The first two terms on the r.h.s describe the change of the superfiuid density 
due to supercurrents between the neighboring layers. In the stationary state 
this is balanced by the conversion of quasi-particle charge into condensate 
charge (the last term in ([T8|)). We assume that in the stationary case this 
conversion can be described by a relaxation process 

dpi 1 ^ 

-g^\con. = --pl (19) 

Using this epression in (|TBp the result ( P^ ) is obtained. 

The relation (p!7|) can also be written in another way which will be useful 
in the following. Using the continuity equation for the total current density 
in the stationary case we find for the quasi-particle current density: 

do'' 

^ Jn—l,n Jn,n+1 ^ q_i^ \convi \'^^ ) 
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then 

'^n = dUZl,n-3Zn+l)/^q (21) 

with cTg = 2elN{Q)d/Tq. 

For the later discussion it will be useful to introduce a relaxation param- 
eter rjn,n+i = o'n,n+i/o'q. In the foUowiug we will assume that nonequilibrium 
effects are small, i.e. a -C 1, r/n.n+i <C 1. In first order in a and rj we may 
then use the approximation j'J^n+i = crn,n+iVn,n+i/d for the calculation of 
and obtain 

— Vn-l,nVn-l,n — Vn,n+lVn,n+l, (22) 

or for a barrier in the resistive state we may replace jn^n+i (^3) by the 
total current density j. 

The relaxation parameter rj can also be written as 77 = Tg/rt, where Tt is a 
tj^ical tunneling time. In the case of weak nonequlibrium, rj <^ 1 considered 
here, the charge imbalance has relaxed before quasi-particles have tunneled 
to the next layer. 

The charge imbalance potential vanishes for a layer between two junc- 
tions in the superconducting state, when the current is carried exclusively 
by Cooper pairs. It vanishes also for a layer between two junctions in the 
resistive state, when the current is carried primarily by quasi-particles. A 
non-zero value of \E'„ can be found, in particular, for a layer between one 
junction in the superconducting and one junction in the resistive state. As 
Tq is proportional to 1/A close to Tc, but jc is proportional to the charge 
imbalance vanishes at Tc. 

These results are derived here for the stationary case, which will be needed 
for the discussion of dc-properties in the following. In the presence of a 
time-dependent chemical potential \E'„(t) an additional term proportional to 
d'^n{t)ldt has to be added in (plOl). The results for the time-dependent case 



are discussed in HIT . 



3 Influence of nonequilibrium eff^ects on the 
IV- curves 

As a first application of the theory outlined above we study the infiuence of 
nonequilibrium effects on the IV-curves, i.e. we calculate the total voltage 
for a stack of junctions, where one or more Josephson junctions are in the 
resistive state in the absence of high-frequency irradiation. Special attention 
is paid to the contact voltage. Here we consider explicitly only one contact. 
The effect of the other contact can be added easily to the final result. 
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3.1 All junctions in the superconducting state 

Let us start the discussion with the case where all junctions are in the super- 
conducting state. Then the voltages are determined by the following set of 
equations: Equation ( [T^ ) is used for the contact with the normal electrode 



and ([TT]) with n > 1 for the junctions in the superconducting state. Adding 
up these equations we obtain for the total voltage which is the contact voltage 
in the superconducting state: 

V = Kent = ^ + ^1. (23) 

For the quasi-particle potentials we have: 

^1 = Vo,ijd/(ro,i = jd/ag, ^„ = for n> 2. (24) 
Then the contact voltage is 

Kent = Kent + SV (25) 

with SV = jd/(Tg and the "bare" contact voltage V^^^^ = jd/ao^i. We see 
that the contact voltage is changed by the charge imbalance induced by the 
quasi-particle current on the first superconducting layer. 

3.2 One or two junctions in the resistive state 

If the first Josephson junction (between layers 1 and 2) is in the resistive 
state, but all the other junctions are in the superconducting state, then the 
contact with the normal electrode is again described by (p^. For the resistive 
junction we use the current relation (15) with n = 1, for the superconducting 



junctions (|TT]) with n > 2. 

Adding up these equations we find for the total voltage: 

V = ^ + ^ + ^2, (26) 
0"o,i cri,2 

where we have used \E'„ = for n > 3. The nonequilibrium effect now comes 
from the second layer, where = Vi,2jd/(Ji^2 = jd/aq. Here and in the 
following we assume that the charge imbalance relaxation rate is the same 
on all superconducting layers. Then the total voltage can also be written as 

V = Kent + — . (27) 

0^1,2 
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The IV-curve of the first branch appears at the usual voltage after subtracting 
the contact voltage. This is different, if the Josephson junction in the resistive 
state is not close to the normal contact. Then one obtains for the total 
voltage: 

V = Kont + -^(1 + 2r7„,„+i) = Kont + + 5V (28) 

with 5V = 2jd/aq. The first branch has a slightly larger voltage for the 
same current density. The voltage shift 6V is due to the charge imbalance 
potentials generated by the quasi-particle currents on the two layers n and 
n + 1 between the resistive junction and the neighboring superconducting 
junctions. 

Of special interest is the case where two junctions are in the resistive state, 
if we compare the total voltage for the case where the resistive junctions 
are neighbors with the situation where they are separated by one ore more 
barriers in the superconducting state. In the latter case the total voltage 
is again larger by 6V = 2jd/aq. This shift is due to the charge imbalance 
potential on the two additional boundaries separating the resistive junctions 
from the junctions in the superconducting state. This result can be easily 
generalized to the case of several resistive junctions. 

In the present theory the shift of the chemical potential described by the 
coefficient a has no influence on the total voltage. This is a consequence of 
the special form of the quasi-particle current density used here (compare with 
([T^ , 0)). This result is different from earlier treatments [0, where for the 
quasi-particle current density a simple ansatz of the form = aVn,n+i/d 

has been used. 



4 Influence on Shapiro steps 

For a Josephson junction in the presence of strong high-frequency irradiation 
the current relation (ITH) is not useful, since on a Shapiro step the current is 
not fixed by the voltage. Instead the Josephson oscillation frequency is locked 
to the external radiation. Then we have to use the generalized Josephson 
relation equation (0) to determine the voltage. Let us assume that the first 
Josephson junction (between layers 1 and 2) is in the resistive state and 
all other junctions are in the superconducting state, then the voltages for 



the different barriers are again determined by (16, 11), but for the resistive 
junction with a first Shapiro step we have to use 

^ = (1 + 2a)Vi,2 - a(\/o,i + ^^2,3) + ^2 - ^i- (29) 
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With ^„ = for n > 3 we obtain for the total voltage: 

jd fiuj ^ , , 

y = — + — + *!. (30) 
cro,i ie 

In this case the quasi-particle potential \E'i is given by 

'^i-d{jl^,-jf;,)/a,. (31) 

As jofi = 3 and jf^^ ~ j at the center of the Shapiro step, vanishes in this 
case. Therefore we find for the total voltage 

V^ = KL + ^, (32) 

where V^^^^ — jd/ao,i is the "bare" contact voltage. If we introduce instead 
the contact voltage which is measured in the superconducting state we have: 

V = Kont + ^ - 5^ (33) 

with SV = jdrio,i/ao,i = jd/aq. 

The apparent shift of the Shapiro step for the first Josephson junction 
is due to the difference of the charge imbalance potential ^'i on the first 
superconducting layer in the pure superconducting state and in the resistive 
state (the effect of the charge imbalance ^2 on the second layer drops out in 
the total voltage). In order to get such an experimental shift of the Shapiro 
step voltage it is essential that the Josephson junction in resonance with 
the external radiation is close to the normal electrode. Otherwise the quasi- 
particle potentials drop out in the total voltage. 

The result can easily be extended to the case where the first m Joseph- 
son junctions are in the resistive state and in resonance with the external 
radiation. In that case we find 

V^V,ont + m^-SV (34) 

with 6V = jdrjo i/cTQ i = jd/cTg. Again it is assumed that the first Josephson 
junction in the resistive state is close to the normal electrode, otherwise 
SV — 0. Note that only the voltage of the Shapiro step on the first branch is 
shifted, the difference between the other Shapiro steps is still at the "right" 
voltage huj/ (2e). 

A similar result is obtained for higher order Shapiro steps on the first 
branch. Then the voltage of the m*'^ Shapiro step is given by 

V = yeont(j) + mj) = KL(i) + m^- (35) 
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Figure 2: Stack with two normal electrodes used to measure the charge 
imbalance relaxation rate. 

As the Shapiro steps appear here at different current values we have to tak;e 
into account the current-dependence of the contact voltage explicitly. Then 
all Shapiro steps appear to be shifted if the contact voltage in the supercon- 
ducting state is subtracted. 

5 4-point measurements 

In the preceding section we have investigated the influence of charge imbal- 
ance on the shift of Shapiro steps. This effect is indirectly due to a change of 
the contact voltage between the normal electrode and the first superconduct- 
ing layer caused by charge imbalance. With help of a 4-point measurement 
it is possible to determine the contact voltage separately. Furthermore a 
direct determination of the charge imbalance and its relaxation rate will be 
possible. 

We consider a stack with two normal electrodes on top (see figure 2). 
Then one contact can be used to inject a current and the other contact to 
measure the voltage in the absence of a current. 



Let us denote the current densities and voltages of the left and right 
electrode (with length /^..i?) by Jl^r and Vl^r. We assume that the currents 
injected through the left and right electrode are different, in particular we 
are interested in the case jR = 0. 

First we consider the case that all inner junctions are in the supercon- 
ducting state. Then the individual voltages are determined by the following 
set of equations: 

^ = (1 + a)V,Y - aVl'f, (36) 

(^0,1 



13 



= (1 + 2a)Vj^f+, - a(Kfif,, + + - *n, for n > 1. (37) 

The injected currents induce a charge imbalance potential on the first 
superconducting layer controlled by the quasi-particle currents. We assume 
that the diffusion length of quasi-particles along the layer is large compared to 
the width of the stack. Then we obtain a uniform charge imbalance potential: 

l.^i'-fh+'-fjR)-. (38) 

I I (7q 

Note that the total charge on the first superconducting layer in general is not 
uniform and the voltages Vl and Vr are different, but the electrochemical 
potential which is controlled by (7n,n+i) is constant along the layers. 
For the voltages at the left and right contact we then obtain: 

Vl,r = ^ + (39) 

(^0,1 

In the case Jr = we find in particular: 

Fl-^k = — , Vr = vMvo,i. (40) 

(^0,1 t 

which allows us to measure the relaxation parameter 770,1 and the "bare" 
contact voltage jid/co,! separately. 

Now we consider the case that the barrier between the first two super- 
conducting layers is in the resistive state. Then a similar calculation leads 
to 

V,,, = ii^ + J± + t, (41) 

0-0,1 0-1,2 

with ^'2 = jdr]i^2/o'i,2 = jd/aq. Now the nonequilibrium effect comes from 
the second superconducting layer. Note that the current density flowing 
through the barrier (1,2) is the total current density j. In the special case 
Jr = a comparison between Vl and Vr allows us to determine the "bare" 
contact voltage of the left electrode. 

Finally let us discuss the results for applied high-frequency radiation. If 
the barrier between the flrst two superconducting layers is in the resistive 
state and is in resonance with the radiation frequency cu the first Shapiro 
step appears at the voltage 

y JM4^r^ + ^^ (42) 

(7o,i Ze 

with ^1 ~ in the center of the Shapiro step. In particular, for = 
the voltage at the right electrode measures the Shapiro step voltage directly 
without the contact voltage. The same is true for the voltage of the m*^ 
Shapiro step on the first branch. 



14 



6 Summary 



In this paper we have discussed the influence of nonequihbrium effects on 
stationary properties of Josephson contacts in layered superconductors. In 
particular, we have studied the influence on the IV-characteristics and on the 
voltage-position of Shapiro steps. We find that from the two basic nonequi- 
hbrium effects, a) shift of the chemical potential of the condensate and b) 
charge imbalance of quasi-particle charges, only the latter has an influence 
on these properties, while the chemical potential shift drops out in the to- 
tal voltage. This result depends crucially on the form of the quasi-particle 
current density and the correct treatment of the contact resistance. Note 
that for time-dependent processes, like collective modes, also the coupling 
between superconductng layers due to the chemical potential shift is relevant 

i- 

In our theory the charge imbalance relaxation plays an important role. 
Microscopically this relaxation is caused by inelastic scattering processes 
within the layer. In the case of d-wave superconductors also elastic scatter- 
ing may contribute. In the present treatment a phenomenological relaxation 
rate has been introduced. 

In our investigation of Shapiro steps we find a shift of the voltage of 
the first Shapiro step on the first resistive branch of the IV-curve if the 
Josephson junction in the resistive state is close to the normal electrode. 
This apparent shift, which is not a violation of the basic Josephson relations, 
can be traced back to a change of the contact voltage between the normal 
electrode and the first superconducting layer due to charge imbalance on the 
first superconducting layer in the presence of a finite quasi-particle current. 
The distance between higher order Shapiro steps on the first branch is also 
changed due to the current dependence of the contact voltage. 

Recently Shapiro step experiments have been performed in the THz regime 



on mesa structures of BSCCO fl^. The voltage-position of the Shapiro step 
has been determined by subtracting the contact voltage in the superconduct- 
ing state in the absence of radiation. In one experiment a 3% downshift of 
the step position on the first branch has been observed, i.e. the ratio be- 
tween charge imbalance lifetime and tunneling time is 77 = Tq/rt = 0.03. In 
other samples also Shapiro steps on higher branches could be detected. Their 
positions showed a more complicated structure which up to now is not fully 
understood. Therefore it is not yet possible to present a detailed comparison 
with experimental data. 

The difficulty with the proper subtraction of the contact voltage can 
be avoided in 4-point measurements. Therefore we have investigated the 
influence of charge imbalance on the Josephson effect for an experimental 
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set-up where the first superconducting layer is in contact with two normal 
electrodes, and only one electrode is used to inject the current. Also in 
this case the contact voltage is influenced by the charge imbalance produced 
by the quasi-particle currents, but the contact voltage can be determined 
separately by comparing the voltages measured at the two electrodes. Taking 
that into account, the Shapiro steps appear again at the right position. 
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